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INTRODUCTION 


This  paper  presents  an  applied  method  for  the  optimal  design  of  passive  vibration  absorbers 
to  reduce  terrain-induced  vibrations  of  tank  cannons.  It  is  anticipated  that  this  will  improve  the 
accuracy  of  the  weapon  by  reducing  variations  in  the  initial  conditions  of  the  gun  barrel  at  the 
commencement  of  firing  dynamics.  The  paper  focuses  on  the  attenuation  of  vibration  in  the  vertical 
plane;  application  to  the  horizontal  modes  follows  in  analogy. 

Flexural  vibration  of  large  caliber  gun  barrels  has  become  more  of  an  issue  due  to  the 
increased  exit  velocity  and  fire-on-the-move  demands  placed  on  the  weapon  system  to  defeat 
reactive  armor-protected  threats  at  extended  ranges.  The  velocity  requirement  results  in  longer 
barrels  that  have  lower  fundamental  modes  of  vibrations  and  that  are  susceptible  to  flexural  vibration 
both  before  and  during  firing  dynamics.  Reducing  the  flexural  state  of  the  barrel  caused  by  terrain- 
induced  vibrations — exclusive  of  the  firing  event — via  a  passive  absorber  enhances  the  overall 
performance  of  the  weapon.  Controlling  vibrations  during  the  firing  cycle  is  more  challenging 
because  of  the  speed  of  the  process;  it  is  too  fast  for  an  absorber  to  shift  or  dissipate  significant 
energy.  Thus,  the  goal  is  to  reduce  the  problem  before  shot-start. 

Applying  vibration  absorbers  to  beams  entails  coupling  a  damped  mass-spring  system  to  the 
beam  at  the  location  of  greatest  vibration  activity.  This  achieves  two  main  benefits.  First,  the 
addition  of  the  absorber  reshapes  the  receptance  of  the  system  by  shifting  the  resonant  modes  and 
zeros.  Reshaping  the  frequency  response  away  from  known  disturbance  frequency  bands  effectively 
rejects  the  disturbance  energy.  Second,  the  absorber  enhances  the  dissipation  of  vibrational  energy. 
This  second  benefit  is  commonly  misunderstood  to  always  be  the  dominant  mechanism  by  which 
the  absorber  reduces  vibrations. 

Actively  controlling  structural  vibrations  via  an  improved  elevation  servo-stabilization 
system  is  currently  being  investigated  by  the  Advanced  Drives  and  Weapon  Stabalization 
Laboratory,  ARDEC,  AMSTA-AR-FSF-R.  However,  this  approach  is  complicated  by  the  non¬ 
minimum  phase  (right-hand-plane  zeros)  behavior  of  the  essentially  cantilevered  beam  system  and 
limitations  of  the  current  hardware,  which  may  require  a  significant  upgrade  for  successful 
implementation.  Incorporating  a  passive  absorber  does  not  preclude  integration  with  active  control 
techniques  and  may  even  enhance  them. 

The  method  presented  in  this  paper  uses  a  Euler-Bemoulli  finite  element  technique  to 
generate  the  second-order  equations  of  motion  of  the  gun  barrel  as  a  non-uniform  beam,  with 
subsequent  conversion  to  the  first-order  state-space  domain.''^  This  model  is  then  transformed  to  the 
frequency  domain  using  the  MATLAB®  software  package.^^’^'  The  design  is  optimized  by  assigning 
a  scalar  cost  function  to  the  frequency  response  function  of  the  modified  barrel;  this  provides  a 
metric  for  minimization  in  the  design  parameter  space  of  the  vibration  absorber.  This  applied 
approach  forms  a  middle  ground  between  mere  numerical  simulation  and  analytic  formulation  of  the 
problem. 
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THE  DYNAMIC  MODEL 


The  Barrel 


The  finite  element  method  was  chosen  to  dynamically  model  the  barrel.  This  application 
uses  the  Euler-Bemoulli  beam  approximation  and  Hermite-cubic  interpolation  functions  to  form  the 
inertial  and  stiffness  matrices  of  the  undamped  second-order  equations  of  motion^'"^  by 
approximating  the  continuous  non-uniform  beam  as  a  finite  number  of  discrete  elements.  Within 
each  element,  the  interpolation  functions  are  used  to  approximate  the  interior  deformation.  At  the 
boundary  between  two  elements  (called  a  node),  continuity  of  lateral  displacement  and  slope  are 
imposed.  When  assembled,  the  resulting  finite  element  model  dynamics — governed  solely  by  the 
nodal  displacements  and  rotations — closely  approximate  the  dynamics  of  the  non-uniform  beam. 


Barrel  Profile,  Coupled  Mass  Distribution,  External  Coupling  Points,  and  Remaining  Node  Locations 


C3 

Pi 


Figure  1.  XM291  tank  gun  geometry  and  node  locations  used  for  finite  element  formulation 


The  XM291  geometry  is  shown  in  Figure  1.  The  plot  depicts  the  inner  and  outer  radii  of  the 
barrel  with  respect  to  axial  position.  The  distribution  of  the  coupled  mass  of  components  (such  as 
the  breech,  cradle,  and  recouperators)  attached  to  the  barrel  is  represented  by  the  thickness  of  the 
coaxial  lines  to  the  rear  of  the  barrel.  This  mass  increases  the  effective  linear  density  of  the  beam. 
(Components  (such  as  the  breech)  whose  center  of  mass  lies  beyond  the  beam  geometry  are  explicitly 
coupled  to  the  end-node  of  the  model  after  finite  element  formulation.)  Impact  testing,  mount- 
constraint,  and  vibration  absorber  locations  are  also  indicated.  The  muzzle-end  sensor  location  lies 
to  the  immediate  left  of  the  absorber.  Nodes  are  explicitly  placed  at  these  five  points  to  facilitate 
external  coupling  with  a  10  element  model  of  the  barrel.  Of  the  remaining  six  nodes,  one  must  be 
placed  at  the  muzzle-end  of  the  beam,  which  leaves  five  to  be  evenly  placed  by  a  meshing  metric.^*^ 
Mount  constraints  are  later  applied  to  eliminate  the  rigid  body  modes  of  the  barrel  and  to  emulate 
the  boundary  conditions  of  the  barrel  as  tested.  Finite  element  modeling  of  the  barrel  results  in  the 
unconstrained  22  x  22  inertial  (M)  and  cross-sectional  stiffness  (K)  matrices  that  govern  the  22  x 
1  generalized  coordinate  (q)  and  force  (f)  vectors.  The  generalized  coordinate  and  force  vector 
elements  correspond  to  the  alternating  lateral  and  rotational  states  of  the  1 1  nodes.  Later  on,  these 
matrices  and  vectors  will  be  modified  and  incorporated  into  the  equation  of  motion,  equation  (2). 
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Constraint/Mounting  of  the  Barrel 


Once  the  dynamics  of  the  distributed  parameter  system  of  the  barrel  are  modeled,  they  must 
be  constrained  by  an  approximation  of  the  gun  mount  support.  This  constraint  is  applied  to  the 
barrel  at  two  locations — the  elevation  mechanism  and  the  trunnions.  For  this  paper,  the  test  structure 
was  grounded  by  affixing  a  gun  mount  support  to  a  pair  of  I-beams,  which  were  bolted  and  epoxied 
to  a  concrete  floor  as  depicted  in  Figure  2.  The  gun  mount  support  directly  constrains  the  trunnion 
bearing;  a  steel  shaft  (not  visible)  is  pinned  between  the  elevation  mechanism  attachment  point  and 
the  gun  mount  support  behind  the  trunnions. 

Both  constraint  forces  are  modeled  as  stiff  springs.  For  this  analysis,  both  constraints  were 
tuned  (using  a  gradient  descent  method — the  MATLAB®  <fmins>  command to  match  the  first  two 
experimental  modal  frequencies  of  the  gun  system.  The  equivalent  spring  at  the  elevation 
mechanism  was  set  to  65  x  10^  N/m,  and  the  trunnion  rate  was  set  to  214  x  10*  N/m.  This  altered 
the  previous  frequency  values'^^  by  less  than  10%  and  is  consistent  with  the  stiffer  effective  mounting 
of  the  trunnions  relative  to  the  elevation  mechanism.  Using  arguments  developed  in  Reference  [1], 
the  constraint  stiffnesses  are  added  to  the  diagonal  elements  of  the  finite  element  stiffness  matrix  that 
correspond  to  the  lateral  generalized  coordinates  of  the  constrained  nodes,  Kj  3  and  Kj  5,  respectively. 

This  parametric  formulation  of  the  constraints  is  limited;  the  real  stmcture  would  require 
more  than  two  values  for  an  accurate  model.  This  approximation  results  in  imprecise  boundary 
conditions  on  the  barrel.  However,  the  model  does  capture  many  of  the  dynamic  effects  of  interest 
and  can  easily  be  modified  to  incorporate  more  advanced  constraint  approximations — including 
servo-control  dynamics.'®^  In  addition,  tuning  the  boundary  conditions  to  match  the  lowest  two 
experimental  modes  of  the  system  better  characterizes  the  vibrations  in  the  frequency  regime  of 
greatest  interest. 


Figure  2.  Grounded  mounting  of  gun  system 


Damping 

A  damping  matrix  (Cd) — see  equation 
(2) — ^that  introduces  a  force  opposite  in 
direction  and  proportional  to  the  velocity  of  the 
deformation  is  constructed  by  the  Rayleigh 
proportional  damping  approximation  (Cd  =  aM 
+  PK).^^’^’*'  The  critical  damping  ratios  for  the 
first  two  modes  of  the  gun  system  were 
computed  using  the  half-power  method^*^  on  the 
measured  frequency  response  (see  Figure  9). 
The  proportional  damping  coefficients  were 
then  assigned  using  equation  (1).^’  ’^  The  results 
are  shown  in  Table  1.  The  temporal  units  of  the 
damping  coefficients  are  equivalent  to  damping 
(force/velocity)  versus  mass  and  damping  versus 
stiffness,  respectively. 
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Table  1.  Measured  frequencies,  critical  damping  ratios, 
and  the  calculated  Rayleigh  damping  coefficients 


(1),  =  1.25x2%  rad/sec 

Cl  =  18% 

a  =  2.57  s'* 

CO2  =  25.4x27t  rad/sec 

C,  =  5% 

P  =  69.5x10-^8 

a 

P 


1  (O. 


1  0)2^ 


co,C. 


(1) 


Coupled  External  Vibration  Absorbers 

With  the  inclusion  of  the  absorber,  a  new  energy-storing  degree  of  freedom  is  added  to  the 
total  system.  This  requires  a  new  generalized  coordinate,  q23,  to  represent  the  deflection  of  the 
absorber  and  its  time  derivatives  from  its  equilibrium  position. '''  This  requirement  increases  the  size 
of  the  system  matrices  (to  be  used  in  equation  (2))  from  22  x  22  to  23  x  23. 

Figure  3  depicts  the  population  of  the  completed  system  mass  and  stiffness  matrices  of  an 
absorber  modified  system.  (The  damping  matrix  formulation  is  analogous  to  the  stiffness  matrix.) 
The  plots  verify  the  cascading  construction  of  4  x  4  elemental  matrices — resulting  in  the  diagonally 
banded  system  matrices  for  the  22  x  22  finite  element  portion.  Furthermore,  coupling  the  absorber 
with  parameters  Mva,  Cva,  and  Kva  to  the  lateral  generalized  coordinate,  q,^,  increases  Kjg  ,9  by  Kva; 
places  Kva;  into  K23  23i  inserts  -  Kva  ^,9  33  and  K23 19;  and  places  Mva  at  M23  23. 


Population  of  K  Population  of  M 
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Figure  3.  Population  of  absorber 
modified  system  matrices 
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Equation  of  Motion 


Once  the  combined  system  matrices  are  formed,  the  resulting  equation  of  motion  is  as 
follows; 


Ma  +  +  Kg. 


(2) 


where: 


M  is  the  23  x  23  mass  matrix.  g 

Cj)  is  the  23  x  23  damping  matrix.  £ 

K  is  the  23  x  23  stiffness  matrix.  • 


is  the  23  X  1  generalized  coordinate  vector, 
is  the  23  x  1  generalized  force  vector, 
denotes  d/dt. 


Conversion  to  Single-Input/Single-Output  First-Order  State-Space 

The  second-order  symmetric  equations  of  motion  of  equation  (2)  are  transformed  to  first- 
order  state-space  using  the  method  presented  in  References  [1]  and  [4].  The  multi-input/multi-output 
state-space  model  is  then  truncated  to  a  single-input/single-output  system  for  frequency  response 
analysis.  The  input,/,,  is  located  at  the  node  at  the  rear  of  the  barrel  where  the  modal  hammer 
impacts  are  to  be  applied  (see  Figure  1).  The  measured  deflection  at  the  node  before  the  absorber, 
q,7,  is  selected  as  the  sole  output  of  the  system. 

Bode  Analysis 

Once  in  the  single-input/single-output  state-space  realization,  the  frequency  response  of  the 
system  may  be  computed  using  the  MATLAB®  <bode>  command. The  frequency  response  of  a 
dynamic  system  indicates  the  steady-state  response,  y(t),  of  the  system  to  a  sinusoidal  input, 
as 


u(t)  =  A  sin{(2nf)t) 
y{t)  =  A:i4  sin((27i:/)r  +  4)) 


where: 


u(t)  is  vertical  input  at  the  breech,/  (N). 
y(t)  is  response  of  the  muzzle,  qj^  (m). 

A  is  amplitude  of  the  input  force  (N). 
k  is  the  receptance  gain  (m/N). 

/  is  frequency  of  the  input  force  (Hz), 
(f)  is  phase  of  the  response  (rad). 

t  is  time  (s). 
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FREQUENCY  DOMAIN  OPTIMIZATION  OF 
THE  ABSORBER 

Bode  Analysis  for  Absorber  Optimization 

The  frequency  response  of  the  unmodified 
gun  system  is  compared  to  the  response  when 
various  parametric  absorber  configurations  are 
coupled.  For  this  paper,  two  absorber  parameters 
are  varied — the  stiffness,  Kva,  and  the  stiffness 
proportional  damping  coefficient,  Pva — while  the 
mass,  Mva,  is  held  constant  at  18  Kg  and  the  mass 
proportional  coefficient,  ttyA,  is  set  to  zero.  The 
gain  function  of  the  modified  barrel  is  then  divided 
by  the  gain  of  the  original  system — resulting  in  a 
relative  frequency  response,  Gr(/),  as  shown  in 
Figure  4. 


Figure  4.  Juxtaposition  of  modified  and  un-  ^  weighting  function,  is  applied  to 

modified  barrels  and  the  resulting  relative  emphasize  known  disturbances.  This  function,  as 

trequencyrespoi.se  5 

was  generated  as  the 
normalized  response  of  a  second-order  system 
with  a  natural  frequency, /„,  of  nine  Hz  and  a  critical  damping  ratio,  C,  of  40%.  This  was  chosen  to 
highlight  a  known  but  poorly  characterized  disturbance,  which  is  introduced  to  the  system  through 
the  tank  suspension  (a  low-pass  filter),  and  to  reduce  emphasis  on  the  frequencies  below  the  first 
mode  of  the  system,  which  are  adequately  handled  by  the  existing  elevation  servo-stabilization.  In 
general,  the  weighting  function  combines  the  dynamics  of  the  system  with  any  a  priori  frequency 
characterization  of  the  disturbance.  (Multiplication  in  the  frequency  domain  is  equivalent  to 
convolution  in  the  time  domain'®’).  Thus,  this  weighting  function  approximates  a  white  force 
spectrum  that  passed  through  a  second-order,  damped  filter. 


Weighting  Function  Frequency  Response 


Frequency  (Hz) 


Figure  5.  Weighting  function  used  to 
emphasize  expected  disturbance  frequency 
content  in  the  optimization 
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A  single  scalar  cost  function,  J,  is  then  quantitatively  computed  as  the  normalized  integral 
of  the  weighted  relative  frequency  response  of  the  absorber  modified  system  across  its  effective 
bandwidth 


J 


{G^if)Gyff))df 


G^df 


(4) 


where: 

is  lower  bandwidth  frequency, 
is  upper  bandwidth  frequency. 

Other  cost  functions  could  be  used  to  incorporate  other  design  considerations  (such  as  worst- 
case  performance)  or  pragmatic  engineering  issues  (such  as  bounds  on  reasonable  damping  levels). 

The  flexibility  provided  by  incorporating  virtually  any  weighting  function  (including  servo- 
control  characterization  and  the  potential  for  multi-mode  receptance  reduction)  motivated  the 
development  of  this  technique.  Analytical  techniques  for  the  optimal  design  of  absorbers  that  are 
based  on  minimum-slope  cross-over  points^’®'  or  a  quadratic  cost  function  '‘^’were  not  used  because 
they  do  not  seamlessly  integrate  with  this  formulation. 

Optimization  Surface 

Figure  6  shows  the  optimization  surface  generated  for  an  18  Kg  vibration  absorber  across  the 
bandwidth  0  to  16  Hz  using  the  scalar  frequency  response  cost  function.  The  plots  clearly  indicate 
that  for  low  levels  of  damping,  Pva  <  0-01  s,  the  absorber  design  is  sensitive  to  parametric  variation. 
This  is  seen  by  the  high  density  of  the  contour  lines  near  the  peaks,  which  indicate  step  gradients. 
The  peaks  are  caused  by  detrimental  interaction  of  the  absorber  with  the  fundamental  frequencies 
of  the  barrel  and  numerical  integration  errors  for  the  undamped  cases. 

The  contour  plot  indicates  that  the  optimal  undamped  design — indicated  by  point  B  at  Kva 
=  7,333  N/m — ^would  effect  disturbance  rejection  with  a  J  value  of  0.99.  If  decoupled  from  the 
barrel,  it  would  have  a  natural  frequency, /„,  of  3.18  Hz.  For  higher  damping  levels,  0.02  s  <  P  va> 
the  absorber  design  appears  to  be  quite  insensitive.  The  optimal  absorber,  indicated  by  point  A  at 
Kva  =  22,667  N/m  and  Pva  =  37.9x10'^  s,  would  effect  disturbance  rejection  by  the  absorber  with 


J  is  scalar  cost  function. 

Gr(/)  is  relative  frequency  response. 
Gw(/)  is  weighting  function. 
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a  J  value  of  0.95.  If  decoupled  from  the  barrel,  it  would  have  a  damped  natural  frequency, /j,  of  4. 17 
Hz  and  a  critical  damping  ratio,  C.  of  67%. 

The  lack  of  sensitivity  to  design  parameters  for  an  optimal  absorber  is  a  desirable  result. 
Engineering  obstacles  with  damping  materials  that  are  subject  to  harsh  variations  in  temperature — as 
would  be  experienced  in  a  weapons  environment — may  be  relaxed  because  of  this  lack  of  sensitivity. 


Surface 


Figure  6.  Optimization  surface  in  the  design  parameter  space  of  an  18Kg  absorber 


Furthermore,  the  negative  slope  of  the  line  of  optimal  Kva's  with  respect  to  Pya  is  in-line  with  the 
temperature-dependent  properties  of  visco-elastic  polymers,  one  candidate  technique  for  achieving 
the  absorber  damping. 

The  frequency  response  function  and  relative  response  of  both  of  these  optimal  absorbers  are 
shown  in  Figures  7  and  8,  respectively.  (A  dB  is  related  to  the  gain  of  equation  (3)  as  201og,o(k/kp), 
with  the  reference  receptance,  k^,  set  to  1  m/N.) 

EXPERIMENTAL  VALIDATION 

Experimental  Frequency  Response  Functions 

Two  absorber  configurations  were  experimentally  evaluated  and  compared  to  the  model.  For 
pragmatic  reasons,  both  absorbers  were  relatively  low  in  damping  compared  to  the  optimal  absorber 
design  indicated  in  Figure  6  because  proof  of  principle  demonstration  could  be  achieved  with  simple 
and  available  hardware. 
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Simulated  Frequency  Response  with  Optimal  Absorbers  to  Vertical  Breech  Impact 


Frequency  (Hz) 


Figure  7.  Simulated  frequency  response  function  of  optimal  absorbers  of  Figure  6 


Simulated  Relative  Response  with  Optimal  Absorbers 


Figure  8.  Simulated  relative  frequency  response, 
Gr(/),  of  optimal  absorber  configurations 


The  tests  were  conducted  using  a 
Hewlett-Packard  3566A  spectrum  analyzer,  an 
instrumented  3-pound  PCB  modal  impact 
hammer,  and  Electro-Mike  PA12D  Vz"  eddy 
probe  displacement  sensors.  Ten  averages 
were  used  to  smooth  spurious  variations  in  the 
measured  frequency  response  function.^*' 

The  three  frequency  response  functions 
measured  are  shown  in  Figure  9,  with  the 
relative  response  depicted  in  Figure  10.  The 
absorber  mass  was  set  using  calibrated 
weights.  The  remaining  absorber  parameters 
were  inferred  from  the  power  spectrum  and 
the  transient  decay  of  each  absorber  coupled  to 
a  rigid  stand.  The  damped  natural  frequency, 
/j,  corresponds  to  the  peak  of  the  power 
spectrum.  The  damping  was  evaluated  using 
the  logarithmic  decrement  method.'®  '®'  The 
resulting  absorber  parameters  are  shown  in 
Table  2. 
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Table  2.  Measured  parameters  of  the  two 
18.4  Kg  absorbers  tested 


Kva  (N/m) 

PvA 

/d(Hz) 

C(%) 

11,900 

0.000243 

3.88 

0.296 

22,000 

0.00230 

5.50 

3.97 

Figure  9  clearly  demonstrates  that  the  first  two  modes  dominate  the  response. .  Both  vertical 
modes  suffer  some  degree  of  cross-talk  with  the  horizontal  mode,  which  has  a  slightly  higher 
frequency  because  of  the  rotational  resistance  of  the  opposed  trunnions,  which  is  not  present  in  the 
vertical  mode.  The  plot  also  reveals  a  lack  of  an  asymptotic  stiffness  line'*’  at  the  low  frequencies 
and  a  discrepancy  in  the  dB  scale  relative  to  that  of  the  model.  The  cause  of  these  two  shortcomings 
is  not  certain  and  is  believed  to  be  a  discrepancy  in  the  analyzer  configuration. 


Measured  Frequency  Response  with  Test  Absorbers  to  Vertical  Breech  Impact 


Figure  9.  Experimental  frequency  response  function  of  unmodified  gun  system  and  two 
absorber  configurations 


Simulation  of  Tested  Absorber  Configurations 

Figures  11  and  12  depict  the  simulated  frequency  response  of  the  three  hardware 
configurations  tested. 

A  comparison  with  Figures  9  and  10  indicates  that  the  poles  and  zeros  of  the  absorber  were 
reasonably  placed  by  the  model.  Discrepancies  in  the  peak  amplitudes  of  the  absorber  response  are 
expected.  Both  absorbers  are  low  in  damping  and  rely  heavily  on  the  dissipation  of  the  barrel  to 
bound  their  response.  These  high  peaks  are  “smoothed”  by  the  finite  fequency  resolution  of  the 
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experimental  results.  The  Bode  analysis  of  the  model  is  not  subject  to  this;  therefore,  the  peak 
response  is  sensitive  to  the  specific  frequencies  evaluated  in  the  vicinity  of  a  pole  or  zero.  It  is 
important  to  note  that  the  cost  function,  J,  is  sensitive  to  variation  in  these  narrow  peaks  if  the 
frequency  vector  used  to  numerically  evaluate  the  integral  of  equation  (4)  is  too  coarse.  The  shift 
to  higher  frequency  poles  than  were  predicted  for  the  stiffer  absorber  may  be  caused  by  the  cross-talk 
between  the  vertical  and  horizontal  modes  of  the  test  structure. 


Measured  Relative  Response  with  Tested  Absorbers 


Figure  10.  Experimental  relative  frequency 
response,  Gr(/),  of  tested  absorber  configurations 


Simulated  Frequency  Response  with  Tested  Absorbers  to  Vertical  Breech  Impact 
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CONCLUSION 


Simulated  Relative  Response  with  Tested  Absorbers 


Applying  a  gun  barrel  vibration 
absorber  to  favorably  reshape  the  frequency 
response  of  a  tank  cannon  was  shown. 

Significant  reductions  in  amplitude — greater 
than  50% — were  demonstrated,  both 
experimentally  and  by  the  finite  element 
model. 

A  general  optimization  procedure  for 
designing  an  absorber  that  is  insensitive  to 
parametric  uncertainty  and  that  incorporates  a 
priori  disturbance  characterization  was 
presented. 

The  model  of  the  coupled  absorber/gun 
system  was  compared  to  experimental  results. 

Tuning  the  unmodified  barrel  model  to  match 
the  experimental  fixture  frequencies  and 
damping  ratios  was  implemented.  This  tuning 
does  not  artificially  match  the  absorber  effects; 

instead  it  provides  a  consistent  baseline  from  which  to  evaluate  the  incorporation  of  the  absorber  into 
the  model.  Comparing  the  simulated  and  experimental  frequency  responses  indicates  that  the 
modified  system  poles  and  zeros  were  properly  placed  by  the  model.  Discrepancies  in  the  narrow 
relative  peaks  of  the  absorber  modified  model  response  were  explained  by  the  smoothing  effect 
inherent  in  the  finite  resolution  experimental  frequency  response,  which  is  not  present  in  the  Bode 
analysis  of  the  model. 


Frequency  (Hz) 

Figure  12.  Simulated  relative  frequency  response, 
G(j(/),  of  tested  absorber  conflgurations 
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